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Abstract 

Sensitive dependence of nonlinear systems on initial conditions or parameters 
can be useful in applications. We propose in this paper that bubbling behavior 
in simple driven symmetrical maps may be used as a working principle of 
sensitive devices. The system is stable when there is no input and displays 
bursting behavior when there is small input. The symmetrical property of the 
bursting pattern is very sensitive to the bias of the noisy inputs, which makes 
the system promising for detecting weak signals among noisy environment. 
PACS number(s): 05.45. +b; 
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A common property of many nonlinear systems is their sensitive dependence on initial 
conditions or parameters. This effect can be useful in applications. For example, the sensitiv- 
ity of a chaotic system can be used to control its state to unstable periodic orbits embedded in 
it [1], in targeting the state of the system to desired points in the state space[2], to control the 
system to follow a desired goal dynamics in order to synchronize with another system [3] or to 
allow a message being encoded in a chaotic series for the purpose of secure communication [4], 
only by small modifications of the parameters or state of the chaotic system. The capability 
of achieving quite different behavior by applying only small perturbations improves greatly 
the flexibility of a system to be used in various applications. 

By definition, sensitivity is referred to as the growth of small perturbations to the system. 
So, naively, sensitivity of nonlinear systems can be used to design sensor devices. Many 
systems possess a period-doubling bifurcation when some parameter is varied. Near the onset 
of a period-doubling bifurcation, any dynamical system can be used to amplify perturbations 
near half the fundamental frequency [5]. One disadvantage associated with the application 
of such parameter sensitivity for sensor device purpose is that the control parameter of the 
system must be located extremely close to the critical value of the bifurcation. 

Recently, Bohme and Schwarz proposed to use two identical chaotic systems to construct 
sensitive devices[6]. In particular, they employed the following symmetrically coupled chaotic 
systems 



named chaotic bridge as a sensor device. Sin represents a constant input to be sensed. The 
coupling gain k is chosen near the threshold kc of synchronization, so that for Sin — 0, the 
coupled systems are in synchronization state, and the output Sout = Ik ~ — 0; while 
for Sin 7^ 0, the symmetry of the chaotic bridge is broken, and it may have a large output 
at some moment. Since the synchronization manifold is transversely stable for = 0, 
there must exist local instabilities in the system in order to obtain amplification of small 
perturbations. As pointed out by the authors in [6], in the neighborhood of the boundary kc 
of synchronization, one can expect the highest sensitivity of the system. 

We would like to highlight the connection of the working principle of the above device 
to the phenomenon of attractor bubbling studied recently [7- 10]. When k is just beyond the 
threshold kc, the synchronization manifold is transversely stable. However, there still exist 
some invariant sets, such as the unstable periodic orbits embedded in the synchronization 
manifold, which are transversely unstable. As a consequence, small perturbations in the sys- 
tems which destroy its synchronization manifold will result in large intermittent bursts from 
the synchronization manifold, no matter how small the perturbations are [8]. This is the origin 
of the sensitivity of the above system. The difficulties of application of the system for sen- 
sor devices lie in practical implementations. Just like additive perturbations, any parameter 
mismatches between the systems can also lead to intermittent bursts. Parameter mismatches 
are inevitable in experiment implementations. This is the reason that intermittent desyn- 
chronization was observed beyond the threshold of synchronization in many experiments of 



X = f{x) - k{x -y) + s, 
y = f{y) + k{x - y) - s, 



(1) 
(2) 
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synchronization between well matched electrical circuits[7-10]. This effect imposes great dif- 
ficulties in the experimental implementation of the above sensor devices, because inevitable 
parameter mismatches lead to large output even for Sj„ = 0. The above devices can work 
only if the two systems are ideally identical, which is extremely difficult to realize. On the 
other hand, small external noise can also result in large bursts when S'j„ = 0, which makes it 
very difficult to tell a signal from noise which is always present in the practical environment. 

To avoid the above difficulties, we propose in the following to use simple driven systems 
as sensor devices. 

Attractor bubbling and on-off intermittency[ll,12] are common behaviors that occur in 
coupled nonlinear systems which possess an invariant manifold. They can be achieved in 
very simple parametrically driven one-dimensional maps[12] 

Vn+l = ZnfiVn), (3) 

where /(O) = 0,df{y)/dy\Q ^ 0, and Zn = aXn > is random or chaotic driving signal with 
density function and a is a parameter. For the purpose of the application of the systems 
as sensitive devices, we require the maps to have odd symmetry, i.e. f{—y) = —f{y)- 
The stability of the invariant manifold y — is governed by the linear equation 

yn+l = Znyn: (4) 

which describes the evolution of small perturbations transverse to the invariant line y = 0. 
Here df{y)/dy\o is absorbed into the parameter a. The transverse Lyapunov exponent A of 
the invariant manifold defined as 

1 ^ 

A = lim — In Zn —< In 2; > 

n=l 

determines the stability of the invariant manifold. The critical point Oc at which A = is 
the onset point of on-off intermittency[12]. For Zn being a uniform random driving signal on 
(0, a), pz = 1/a, and <lnz >= Ina — 1 = gives Gc = e. Just above the onset point, a > Oc, 
the random driven system displays universal features of on-off intermittency behavior, which 
are unaffected by the form of the confining nonlinearity[12]. The nonlinearity of the system 
serves to bound or reject the dynamics back towards small values of y after bursts. For a 
below flc, the invariant manifold y = is stable, but the stability is quite weak if a is near 
ttc- Attractor bubbling occurs in the system when there are inputs of perturbations such as 
noise. 

For a sensitive device, it should be stable when there is no input, and is expected to pro- 
duce large outputs when there are small perturbations. In this paper, we are not considering 
the on-off intermittency behavior with a > a^. For the purpose of sensor purpose, we employ 
the bubbling behavior with a < Gc- The sensor system reads 

Vn+l = Znfiyn) + Sin, (6) 

and Sout — y- When there is no input, i.e. Sj„ = 0, y = is a stable solution, and the 
output Sout — 0. Since the critical parameter Oc and the evolution of small perturbations 
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(5) 



are independent of the form of the nonhnearity, one can choose a map which is simple for 
implementation. For example, we employ a piecewise linear map 



f{y) = { y, 



C2 



CI 



(-ci - C2 - y), y< -ci, 

\y\ < ci, 

(ci + c2-y), y > ci, 



(7) 



where Ci and C2 are two parameters, as shown in Fig. 1. 

If there is no input, Sj„ = 0, staring from a small initial condition i/q — p, we have 
Z/n — ' ' ' ZqP = knP- The average order of the factor kn, which can be defined as 

< lnA;„ >, decreases linearly with n as 



Now, suppose there is a small constant input Si„ = p, starting from yo — 0, the evolution of 
the output reads 



Vn = {Zn-1 ■■■Zo-\ h Zn-2Zn-l + + 1)P = Kp, 

if max(|y„|) = |p|(l+aH a") = - "^i' n<nc^ int[ln(ci(a-l)/|p| + l)/ lna]-l 

, where int[a;] is the interger part of a real number x. The average order of kn, < lnA;„ >, 
cannot be obtained analytically as that for Sin = in Eq. 8. A numerical estimation of 
< lnA;„ > is carried out with 10^ samples of kn. Unlike the case Sj„ = 0, it is an increasing 
function of n, as shown in Fig. 2. 

When n > ric, Vn has nonvanishing probability to exceed t/n > ci. Following Bohme and 
Schwarz, a measure s for the sensitivity of the system to constant input can be defined as 



Since max(A;„) = " ^ increases exponentially with n, an infinitely small input value can 
create a finite output. The largest output max|?/„| = aci is due to the confinement of the 

n 

nonlinearity of the map. The sensitivity s of the system thus goes to infinity and is referred 
to as supersensitivity[6]. 

An important question concerning the system is the time needed for small input to produce 
a large output. We examine the time for a small input p to produce for the first time an 
output Sout ^ Ci. The distribution P{N) of N is the following probability 



< lnA;„ >= n < ln2; >= n(lna — 1). 



(8) 



yi = (^0 + i)p = kip, 

y2 = {ziZo + Zi + l)p = k2P, 



(9) 



N-l 



P{N) = Prob( f]y^<c^f]yN> Ci). 



(11) 



i=l 
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N 

By defining the event E'at = f) Hi < Ci and the corresponding probabihty Ajv = Prob(£^7v), 

i=l 

it follows that 

P(7V)=Ajv-Aiv-i, 

which is a function of both a and p. In principle, A at can be evaluated by the joint density 
^{K) of ki, k2, - ■ ■ , k^, namely, 

^ p{zo)p{z,)---p{z^.,) ^ P(fci-I)p(^)---P(fe7) 

^ ' \J\ k,k2---kN-l 

where J is the Jacobian of the transformation between /cj's and z^'s defined in Eq. (9). 
Specifically, one has 

kN-l.dkN-i ki+i-1 dki 
An = dkN p{-, )-, / p( )-. 

J(k2-l)/a ki ki 

where b — Ci/p ii i > ric and b = if i < ric- However, to the best of our knowledge, a 

closed-form solution of the above integral for any N is not available. 

In the following, we are going to carry out some simulations. We specify ci = 1 and 
C2 = 2 in these simulations. In Fig. 3, as an example, the output sequence is shown for 
Sin = 1-0 X 10^'' at a = 2.6. The dashed lines indicate the switch on and off of the constant 
input. The output in the presence of input is a intermittent process, similar to the result of 
the chaotic bridge in [6] . 

In the next simulation, we estimate P{N) for different values of a and p, as shown in 
Fig. 4. It is seen that the distributions peak at rather small values, and after the peak, 
they decrease exponentially. The average time < > for first putting out ijn > ci is also 
evaluated as a function of a and p in Fig. 5(a) and (b), respectively. So, on average, the 
closer the a to Uc and the larger the input p, the quicker the system reaches a large output. 

The simulation results show that this system may be used as a detector for weak signal. 
However, the above discussion is only valid in a noise-free environment. In the practical 
application of the system as a detector, external noise is unavoidable. The system now reads 

Vn+l = Znf{yn) + Sin + Cn, (12) 

where e„ denotes external noise. It is plausible to assume that e„ has vanishing mean value 
and a Gaussian distribution aN{0,l), with a standard deviation a. The behavior of the 
system in the presence of noise is quite different from the noise-free case, because bubbling 
occurs even without a signal Sj„. Very small external noise can also lead to large output of 
the system, as illustrated in Fig. 6 (a) with a = 2.6 and a = 1 x 10^^. In this system with 
a < ttc, bursting behavior always means that there are some inputs to the system, and the 
largest output will be independent of the inputs. The problem now becomes whether we can 
distinguish that the input is a meanful signal or just noise, and moreover whether we can 
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detect any significant signal among the white noise environment. As will be shown in the 
following, this system is quite promising for this task, because the symmetrical property of 
the bursting behavior is very sensitive to the bias of the inputs. 

This sensitivity is due to the odd symmetry of the map. An inspection of Fig. 6(a) reveals 
that the number of large bursts to positive and negative values is quite symmetric for inputs 
of white noise. When a small positive constant input Sm = P = 0.3 x 10~^ is present along 
with the noise, the symmetry is clearly broken, as seen in Fig. 6(b). Note that the constant 
input p is much smaller than the noise level a in this example. This result indicates that the 
symmetry of burst is quite sensitive to the bias of the total inputs Sin + e„ of the system. 
To characterize the symmetry breaking property quantitatively, we introduce the degree of 
asymmetry D^sy as 

= (13) 

where A^_|_ (A^-) is the number of large burst {\y\ > 1) to positive (negative) values during 
a period of observation time T. Because of the symmetry of the map and that of the white 
noise, one can expect that Dasy ~ for white noise inputs. For positive (negative) constant 
inputs in the noise-free case, it is clear that Dasy — 1(~1)- 

A comparison between Fig. 6(a) and (b) also shows that large bursts occur more fre- 
quently when constant input is present with the noise. We define the bursting frequency F 
as 

F ^ (14) 

We expect that F increases with larger constant input p and larger noise level a. 

A measure of the significance of a constant signal among the noise can be the "signal to 
noise ratio" R = p/a. In the following, simulations are carried out to examine the dependence 
of Dasy and F on R for different noise level a and system parameter a. The results shown in 
Fig. 7 are obtained with T = 2 x 10^. The result of Dasy is very interesting: as a function 
of R, Dasy is independent of the noise level a and is not sensitive to parameter a. Whether 
a constant signal embedded in the noise environment can be detected depends only on its 
significance with respect to the noise level. However, a and a have effect on the bursting 
frequency F, as seen from Fig. 7(b). A combination of Fig. 7(a) and (b) makes it possible to 
determine the noise level and R, and thus to estimate the amplitude of the constant signal 
in a observation. 

To get a good estimation of Dasy and F, T should be large enough. In practice, one may 
not expect a signal that stays constant for so long a time. The term constant signal is a 
concept relative to the time scale of the detector, and the time scale can be controlled in the 
implementation. In numerical experiment of this dimensionless system, we can simulate a 
shorter signal (or a "slower" detector) with smaller T, e.g., T = 2000. In this case, Dasy has 
large fluctuations, as shown in Fig. 8(a) where Dasy of 20 realizations of the driving z„ and 
white noise e„ of the system are plotted for each R value. When R is getting larger, more 
points coincide at Dasy = 1- A good way to examine the fluctuation behavior is to construct 
a histogram of Dasy, as shown in Fig. 8(b) for a — 2.6, a — 10~^. The results show that 
even for quite short signal and low it! value, Dasy has very high probability near Ds — 1. An 
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implication of the results is that several detectors can be used at the same time to detect 
and confirm a short and weak signal embedded in white noise. 

As an example of a little more realistic input signal, we present the response of the system 
to the noisy input y4sin(0.003n) + e„, where the noise level is cr = 1 x 10~^ and the amplitude 
of the sine wave is ^4 = 0.3 x 10~^. Both the total input and the output of the system 
are displayed in Fig. 9. The bursting feature reflects the weak wave among the noise quite 
clearly. 

In summary, we have shown that a kind of very simple driven symmetrical maps below 
the onset point of on-off intermittency have two distinguishing features of (i) being stable at 
the invariant state y — and (ii) being sensitive to small input. In practice, the environment 
cannot be noise-free, and the systems exhibit bubbling behavior in the presence of noise. 
Another interesting and useful property of the systems is that the bursting pattern is sym- 
metrical for white noise input, and the symmetry is broken when there is signal among the 
noise environment. The significance of the signal is manifested by the degree of asymmetry in 
the bursting pattern. These features make them promising candidates for designing sensitive 
devices. 

Although our study is based on numerical simulations of a mapping model, it should be 
noted that system response to small inputs is governed by its linearized equation, and the 
nonlinearity only serves to keep the system bounded. Many long time properties shown above 
thus are universal in a class of driven systems possessing odd symmetry. The following can 
be advantages for such systems when considering applications in sensitive devices: 

1) The sensitivity is maintained in a large range of parameter below the critical point. 
This avoids the difficulty of locating parameter in a very small neighborhood of a bifurcation 
point in a period-doubling system. 

2) The sensitivity of the system is, in principle, infinite. In a noise-free environment, an 
infinitesimal input signal can produce a finite output. When noise is present, weak signal 
can also be manifested by the asymmetry in the bursting pattern. The degree of asymmetry 
depends on the significance of the signal with respect to the noise. The sensitive behavior 
is universal for different forms of nonlinearity of the systems, as well as for different form of 
driving signals. This is very useful because one can thus choose a system that is simple and 
easy to implement in practice. 

It could be meaningful to consider implementation of such simple systems and explore 
their application in small signal detection. Since in applications, pivotal role is played by 
the symmetry properties of the system, one should take care to maintain such properties. 
In order to avoid perturbations which may make the system appreciably asymmetric, one 
should avoid using different parameters for the two symmetrical parts of the systems. Also, 
one should note that it takes longer for the system to produce static output states for lower 
level of inputs. There seems to be a frequency cutoff associated with the input levels and the 
relaxation time of the systems. Above the cutoff, the small signal in the noise can no longer 
be manifested by clear asymmetry in the bursting pattern. Such limits should be taken into 
consideration in applications. 
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Figure Captions 



Fig. 1 The piecewise linear map f{y). Here ci = 1 and C2 = 2. 

Fig. 2 < lnA;„ >, the average order of as a function of n. It decreases hnearly for the 
case without input (plot a) and increases for the case with input (plot b). 

Fig. 3 An illustration of the output process of the sensitive system at a = 2.6. The constant 
input is p = 10~^, and is switched on and off alternately for every 500 iterations, as 
shown by the dashed lines. 

Fig. 4 Numerically evaluated distribution of A^. The parameters of the plots are: (a) 
a = 2.5, p = 10-3, (b) a = 2.5, p = 10"^ and (c) a = 2.6, p = 10"^ 

Fig. 5 (a) Average value of A" as a function of a with p — 10"^, and (b) as a function of p 
with a — 2.6. 

Fig. 6 (a) Output of the system when only white noise is present as input, (b) Output of the 
system when constant signal p = 0.3a is present along with the noise. The parameters 
are a = 2.6, cr = 1 x 10 Noting the change of symmetrical property of the bursting 
pattern of the system. 

Fig. 7 (a) Degree of asymmetry Dasy as a function of "signal to noise ratio" R for different 
a and a. (b) Bursting frequency F as a function of R for different a and a. (b) shares 
the same legends of (a). The results are obtained from observation during a period of 
time T = 2 X 10^ 

Fig. 8 (a) Dasy of short time observation, T = 2000. (b) Normalized histograms of Dasy for 
different R value. The histograms are constructed with 50000 observations for each R 
value. The parameters are a = 2.6, cr = 1 x 10~^. 

Fig. 9 (a) A weak sine wave embedded in the noise, (b) The response of the detector to the 
inputs of (a). 
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